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We study three-electron energy spectra in Si/SiGe single and vertically coupled double quantum 
dots where all the relevant effects, such as, the Zeeman splitting, spin-orbit coupling, valley coupling 
and electron-electron Coulomb interaction are explicitly included. In the absence of magnetic field, 
our results in single quantum dots agree well with the experiment by Borselli et al. [Appl. Phys. 
Lett. 98, 123118 (2011)]. We identify the spin and valley configurations of the ground state in the 
experimental cases and give a complete phase-diagram-like picture of the ground state configuration 
with respect to the dot size and valley splitting. We also explicitly investigate the three-electron 
energy spectra of the pure and mixed valley configurations with magnetic fields in both Faraday and 
Voigt configurations. We find that the ground state can be switched between doublet and quartet 
by tuning the magnetic field and/or dot size. The three-electron energy spectra present many 
anticrossing points between different spin states due to the spin-orbit coupling, which are expected 
to benefit the spin manipulation. We show that the negligibly small intervalley Coulomb interaction 
can result in magnetic-field independent quartet-doublet degeneracy in the three-electron energy 
spectrum of the mixed valley configuration. Furthermore, we study the barrier-width and barrier- 
height dependences in vertically coupled double quantum dots with both pure and mixed valley 
configurations. Similar to the single quantum dot case, anticrossing behavior and quartet-doublet 
degeneracy are observed. 



PACS numbers: 73.21. La, 73.22.-f 61.72.uf 71.70.Ej 

I. INTRODUCTION 

Silicon quantum dots (QDs) are proposed to be promi- 
nent candidates for spin qubits owing to the long spin- 
decoherence timefi^— which is of great importance for 
coherent manipulation, information storage and quan- 
tum error correctionJ^^— The long decohercncc time 
in Si-based devices results from the weak hyperfine 
interaction;^ small spin-orbit coupling (SOC)^^"— and 
weak electron-phonon interaction.^^ As an indirect-gap 
semiconductor, the conduction band of the bulk Si has 
six degenerate minima or valleys. This six-fold degen- 
eracy can be lifted by strain or quantum confinement, 
e.g., it separates into a four- fold degeneracy and a two- 
fold one in [001] quantum wells. The two-fold degenerate 
valleys with low energy can be further lifted by a valley 
splitting due to the interface scattering i^^'^^ The exis- 
tence of the valley degree of freedom makes the Si-based 
qubits more attractive.^^i^ Moreover, the mature micro- 
fabrication technology of the classical Si-based electronics 
is proposed to benefit the realization of Si spin qubitsjii^ 

Recently, Si quantum dots (QDs) have been 
widely investigated both experimentally and 
theoretically j^ '^^'^^" — In the experiments. Si metal- 
on-semiconductor and Si/SiGe QDs with a tunable 
electron filling number from zero have been fabricated, 
where the valley splitting, few-electron energy spectrum 
and spin relaxation time have been measured.—""— 
The theoretical works mainly focus on the one-electron 
Zeeman sublevels or the singlet-triplet states of two 



electrons i^""—i^ For example, Culcer et a/.^'^'^-'^ analyzed 
the initialization and manipulation of one-electron and 
two-electron qubits in lateral coupled double quantum 
dots (DQDs) by utilizing the valley degree of freedom. 
Raith et al^ calculated the energy spectrum and the 
spin relaxation time in one-electron QDs. By explicitly 
including the electron-electron Coulomb interaction, 
Wang et al. obtained the two-electron spectrum from 
exact-diagonalization method and calculated the singlet- 
triplet relaxation time in both single^ and lateral 
coupled doublc^° QDs. As pointed out by Barnes et 
al.,^^ the robustness of the quantum states against 
charge impurity and noise in QDs can be improved 
by increasing the number of the electrons, which 
reveals the necessity of the theoretical investigation 
on multi-electron spin qubits. The energy spectra in 
multi-electron GaAs QDs have been explicitly calculated 
to identify the specific spin configuration of each state, 
where only one valley is relevant i^""— However, to 
the best of our knowledge, there is no report on the 
explicit calculation in silicon QDs with three or more 
electrons due to the complication of the calculation. 
Alternatively, Hada et al.^^ neglected the correlation 
effect and calculated the three-electron energy spectrum 
within single configuration approximation in single 
QDs. However, the correlation effect has been shown to 
present significant infiuence on the energy spectrum via 
strong Coulomb interaction in Si QDs^ i'^^'^" Therefore, 
in order to obtain the accurate convergent spectrum, 
the exact-diagonalization method with sufficient basis 
functions is required. The goal of the present work is 
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to analyze the energy spectrum in the three-electron 
Si/SiGe QDs based on the exact-diagonalization method. 

In this work, we calculate the three-electron energy 
spectrum in both single and vertically coupled double 
QDs with the Zeeman splitting, SOC, valley coupling 
and electron-electron Coulomb interaction explicitly in- 
cluded. In the single dot case, we first calculate the 
ground state energy in the absence of the magnetic field, 
where good agreement with the experimental data is 
achieved. Our calculation also uncovers the valley and 
spin configurations of the ground state in experiment. 
We present a complete phase-diagram-like picture to de- 
scribe the spin and valley configurations of the ground 
state. We find that the ground state is of pure (mixed) 
valley configuration at large (small) valley splitting. The 
spin configuration of the ground state can be controlled 
through dot size for the pure valley configuration, while 
that in the mixed valley configuration is always doublet. 
The magnetic field dependence of the three-electron en- 
ergy spectrum in each "phase" is investigated. We take 
into account both orbital effect and Zeeman splitting for 
the perpendicular magnetic field and only Zeeman split- 
ting for the parallel magnetic field owing to the strong 
confinement along the growth direction. We find that 
the spin configuration of the ground state can also be 
switched by magnetic field. For the mixed valley con- 
figuration, we find interesting quartet-doublet degener- 
acy, which results from the negligibly small intervalley 
interaction. In the DQD case, the barrier-width and 
barrier-height dependences of the energy spectrum with 
both pure and mixed valley configurations are discussed. 
Moreover, we show many anticrossing points resulting 
from the SOCs in all cases. 

This paper is organized as follows. In Sec. II, we set up 
our model and formalism. In Sec. Ill, we show our results 
of the three-electron energy spectrum in single QDs and 
vertically coupled DQDs from the exact-diagonalization 
method. We investigate the perpendicular and parallel 
magnetic-field dependences in single QDs and the barrier- 
width and barrier-height dependences in vertically cou- 
pled DQDs. The comparison with experiment in single 
QD case is also given in this section. Finally, we summa- 
rize in Sec. IV. 



II. MODEL AND FORMALISM 

We set up our model in a double quantum well along 
[001]-direction. The confinement along this direction is 
described by^i^ 

r ^0, 1^1 < f, 

V.{z) = I 0, f < |z| < f +d, (1) 
I oo, otherwise, 

with the inter- well barrier height denoting as Vq. Here, d 
and a represent the width of each well and that of the bar- 
rier, respectively. The lateral confinement is chosen to be 



a parabolic potential Vc{x,y) = \mtU}Q{x^ + 2/^) with rrit 
representing the in-plane effective mass and loq being the 
confining potenti al frequency . "'^^^^"^ The effective diameter 
do is given by hm j (rritUJo). The total confinement po- 
tential then can be written as V^(r) = Vc{x,y) + Vz{z). 
For an infinitesimal barrier width (a ^ 0), our model 
reduces to the single dot case. 

The external magnetic field with perpendicular and 
parallel components is expressed by B = B±z, + S||X. 
The single-electron Hamiltonian reads^ 

p2 ^ p2 p2 

Ho = ^ + + ^(i") +Hz + + i?v, (2) 

Zrrit Zniz 

where niz denotes the effective mass along the growth 
direction and P — —ih'V + (e/c)A with the vector 
potential A = {—yB±,xB±,2yB\\)/2. In our calcu- 
lation, the orbital effect of the parallel component of 
the magnetic field is neglected due to the strong con- 
finement along the z-direction and the vector poten- 
tial in the mechanical momentum then reduces to A = 
{—yB±,xB±,0)/2. In Eq. the Zeeman splitting is 
given by Hz = ^5Ms(S_l(t^ + B\\aa:) with g being the 
Lande factor. The SOCs, including the Rashba term^ 
due to the structure inversion asymmetry (SIA) and the 
interface-inversion asymmetry (HA) term,— ^— are ex- 
pressed by 

Hso = aQ{Pxay - Pycr^) + bo{-Px<yx + Py<yy), (3) 

where ao (6o) represents the coupling coefficient of the 
Rashba (IIA) term. in Eq. ^ describes the coupling 
between the two low-energy valleys lying at ±(A;si) along 
the z-axis with (/csi) = 0.85(27r/asi).— Here, asi = 5.43 A 
stands for the lattice constant of silicon. In this work, we 
use "z" ("z") to denote the valley lying at (fcsi) (— (fcsi)) 
for convenience. 

In order to build up a complete set of single-electron 
basis functions, we define i/p = iJo + with i/o = 

+ ^ — I- V{v). Due to the separation of the lat- 
eral and vertical confinement, one can easily solve the 
Schrodinger equation of Hq. From the lateral part, one 
obtains the eigenvalues^"— 

Eni^hQ.{2n+\l\ + \) + hluJB, (4) 

where SI = -^/wq + cj^ and ujb — gBj\_/ (2mt). The eigen- 
functions read^"— 

K^iir.e) = iV„,z(«0'"e-("''^'/'e'"'il!l[(ar)2], (5) 

with Nn,i = {a^n\/['K{n+ and a = ^mtVt/h. 

in is the generalized Laguerre polynomial. Here, 
n — 0,1,2,... is the radial quantum number and I ~ 
0, ±1,±2,... is the azimuthal angular momentum quan- 
tum number. For the vertical part, we include the lowest 
two subbands, one with even parity (denoted by subscript 
riz = 0) and the other with odd parity (n^ = 1). The 
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corresponding wave functions can be expressed as^ 

Cosin[fco(2- f -d)], f < ^ < f + d, 

Aocosh(/3o2:), \z\ < §, 

Co sin[fco(-2 - f - d)], -§ - d < z < -f , 

0, otherwise, 



(6) 



and 



Cisin[fci(z- f -d)], f<z<f+d, 
Aisinh(/3iz), |z| < f, 

- S c'^sin[fci(^ + f + d)], -f-d<z<-f, 



0, 



otherwise, 



(7) 



with fc„^ 



y/2jn^E„Jh? and /3„^ 



\/2mz(Vb — En^)/^^. En__ is the eigenvalue of the 
Uj-th subband. For the single dot case, only the 
lowest subband [uz = 0) is relevant. With the knowl- 
edge of the eigenfunctions of Hq, one expresses the 
single-electron basis functions in different valleys as 
'/'ntal = ■f^"!Cn,e±''=siZy^^_^(j.)^ ^■^^]^ Mz(z)(r) being the 
lattice-periodic Bloch functions. 

Then we introduce the valley coupling accord- 
ing to Ref. m, where the relevant components are 
given byl^S {<l^l,,n.\H.WniJ = {€'l'n^H.\<t>llJ = 

^l,,njn,n-5l,l' and (</>^' i'« J^v | </>^i„ J 

iK'Vn'\HAVnln,) = ,n'Jn.n' Sl,l' . The CXprCSSiouS of 

„^ and A^^ are given in Appendix [K\ Since the 
valley coupling between different subbands {uz ^ n'^) in 
our case is much smaller than the intersubband energy 
difference (< 2 %), we neglect the intersubband valley 
coupling. The eigenstates of Hq then can be written 
as (l>nin, = -^{Vnin, ± CinJ ^nd the Corresponding 
eigen-energies are £'^;„^ = E^i + En, + E^^. with 
= ^n,,n, ± I^L.nJ representing the energy from 
the valley degree of freedom and A£'^^ — 2|A,\ 
being the valley splitting. 

For the three-electron case, the total Hamiltonian can 
be expressed as 



Htot — H} + + + H} 



(8) 



The superscripts "1" , "2" and "3" on the right hand side 
of the equation label the relevant electrons, e.g., rep- 
resents the single-electron Hamiltonian of the i-th elec- 
tron given by Eq. ([2]) and H^-' — 4^gpK'|r -r | stands for 
the Coulomb interaction between i-th and j-th electrons. 
Here, k is the relative static dielectric constant. 

By using the single-electron functions {(pninA '-'^ 
{InlnzHv}} (riv — i denotes the valley eigenfunction) , 
we construct the three-electron basis functions in the 
form of either doublet {Stot = ^, denoted as \Dgl^)) 



or quartet (S't, 



denoted as IQg"^)) with Clebsch- 



Gordan coefficients. Subscript S* stands for the spin 



magnetic quantum number. Superscript S denotes four 
valley configurations of three-electron basis functions, i.e. 
S = — (+) for the case three electrons in "— " ("-f") val- 
ley and S = m (in) for two electrons in "— "("-|-") valley 
and one electron in "-|-"("— ") valley. The details of the 
three-electron basis functions are given in Appendix [B] 

One calculates the matrix elements of Eq. ([5]) under 
the three-electron basis functions. The details of the 
Coulomb interaction are given in Appendix[Cj By ne- 
glecting the small coupling between basis functions of 
different valley configurations,ii^ the complete basis func- 
tions can be divided into four individual sets according 
to valley configurations (S = — , m, m and We diag- 
onalize the Hamiltonian of each subspace and define an 
eigenstate as doublet {D^gt') (quartet IQ^"^)) if its am- 
plitude of doublet (quartet) components is greater than 
50 %. In the following, we still use the notations li^g""*) 

and IQg"'') to describe the spin properties of the eigen- 
states without any confusion. 



III. NUMERICAL RESULTS 

In our calculation, the effective mass mt — 0.19mo 
and niz = 0.98mo with mo representing the free elec- 
tron mass4^ The strengths of the SOCs between the 
states with the same valley index "±" are taken as 
ao = ±6.06 m/s and bo = ±30.31 m/s (Ref.H). The 
Lande factor g — 2 (Ref.l47l) and the relative static di- 
electric constant k = 11.9 (Ref. 39). In single QDs, we 
take 1430 quartets and 3330 doublets to guarantee the 
convergence of the energy spectrum, while in vertically 
coupled DQDs we take 7316 quartets and 16008 doublets 
correspondingly. 



A. Single quantum dots 

1. Comparison with experiment 

Recently, Borselli et al. investigated four Si/SiGe sin- 
gle QDs, which are labeled as Sil, Si2, Si3 and Si4 sepa- 
rately, in the absence of the magnetic fieldfSi They mea- 
sured the tuning voltages for the injection of an addi- 
tional electron into the QDs, i.e., AVq° = y^^<=^(^<=+i) _ 
Vq^" ^"f^^" ^ where ]/^=^(''^<=+i) represents the gate volt- 
age where the electron number changes between N^. and 
A'c+l. Then the addition energy of an incoming elec- 
tron was determined by Afi^^ = hn^+i ~ — c^AVq", 
where fi^^ — E^°^ ~ ^Tot^^ with E^°^. and a being the to- 
tal energy of the A'c-electron ground state and the energy- 
voltage conversion factor. The experimental data of A/ii 
and Afj,2 of the four samples (noted with the superscript 
*) are listed in TableU) 

For a theoretical study, one can obtain the total energy 
of multi-electron ground states and calculate the addition 
energy. However, the well widths and effective diameters 
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Sil 


Si2 


Si3 


Si4 


AEq (meV) 








0.27 


0.12 


d (nm) 


3.945 


3.945 


3.945 


3.926 


Afil (meV) 


4.520 


3.800 


3.916 


4.680 


Afj,\ (meV) 


4.452 


3.828 


4.008 


4.671 


Ri (%) 


1.5 


0.7 


2.3 


0.2 


Aiit (meV) 


3.226 


2.863 


3.146 


3.594 


A/i2 (meV) 


3.318 


2.824 


3.061 


3.613 


R-y (%) 


2.9 


1.4 


2.7 


0.5 


do (nm) 


31.5 


35.3 


34.1 


30.3 




m/fn 


m/m 




m 


"Stot 


1/2 


1/2 


3/2 


1/2 


Deg. 


4 


4 


4 


2 



TABLE I: The comparison between experimental data-'^ and 
our results in four devices, Sil-Si4. AEq and d are the vaUey 
sphtting and half-well width used in the calculation. A/i* 
represent the experimental addition energy and A^i stand for 
the theory results. Ri stand for the relative error between Afi* 
and Afii. The effective diameter do is determined from our 
calculation.-^ The valley index H denotes valley configuration 
of three-electron ground state in each sample. Stot and "Deg." 
represent the total spin and degeneracy factor of the three- 
electron ground state in each sample, respectively. 
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required for quantitative calculation are unavailable in 
Ref.'sT. Therefore, we first derive the well widths from 
Eq. (|A2I) with the reference value suggested by the ex- 
perimental work (d ~ 4 nm)4^ For Si3 and Si4, we use 
the valley splittings AEq from the experimenti^ Since 
the valley splittings of Sil and Si2 are too small to be 
measured in experiment;^ we take them to be zero in 
our calculation. Then, we calculate the ground state en- 
ergies of one-, two- and three-electron cases by treating 
do as a parameter and determine its value by the least 
square method of A/ii and A/i2.^^ In the previous sec- 
tion, we only introduce the frame of the three-electron 
case. For the one-electron case, one needs to diagonalize 
the single-electron Hamiltonian as given in Eq. © , 
while for the two-electron case we follow the frame in 
Ref.^. 

Our results of addition energies as well as the effec- 
tive diameters of QDs are listed in Table H] Good agree- 
ment with experimental data can be observed for all de- 
vices with the largest relative error less than 3 %, which 
confirms the validity and accuracy of our model. In 
fact, one can also estimate the effective diameter of the 
QDs solely from the fitting of A/ii with the experimen- 
tal value Afil. With the effective diameter obtained in 
this way, we recalculate Afi2 and find that the value also 
agrees well with A/ij (within 4 %). We should point out 
that the Coulomb interaction here is very important in 
these quantum systems. Without explicitly including the 
Coulomb interaction, the addition energy is mainly from 
the orbital energy and the theoretical results become far 



FIG. 1: (Color online) (a) Valley and spin configurations of 
the ground state in the absence of magnetic field. The curves 
represent the effective diameter dependence of the valley split- 
ting energy for the emergence of the degeneracy between the 
lowest levels with "— " and "m" valley configurations. For 
each half-well width d, the ground state is of "— " ("m") con- 
figuration above (below) the corresponding curve. The "— "- 
configuration regime is separated into two parts by the ver- 
tical dotted line, according to the total spin Stot. The valley 
configuration and the total spin are labeled as (H, Stot) in 
each regime. The crosses represent the parameters in Sil- 
Si4 in Ref.[3l|. (b) The orbital energy difference and the in- 
verse number of Coulomb energy difference between the lowest 
quartet and doublet states with "— " valley configuration vs. 
effective diameter for d = 2 nm in the absence of the magnetic 
field. 

away from the experimental data. 

2. Ground state configuration 

In Table IH the valley configuration and total spin of 
the three-electron ground states are also listed. As shown 
in the table, the ground state in Si3 is mainly com- 
prised of the quartets with "— " valley configuration due 
to the large valley splitting. This level is four-fold de- 
generate. However, the valley splitting in Si4 is small 
and the ground state is two-fold degenerate doublet of 
"m" -configuration. For Sil and Si2, the states with "fh" 
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valley configuration are degenerate with those with "m"- 
configuration. Therefore, the ground states are four-fold 
degenerate. Since the states with different valley con- 
figurations are decoupled as mentioned abovCf^ we ex- 
plicitly investigate the three-electron spectrum for each 
valley configuration individually in the following. The 
relative position between the spectra of different valley 
configurations is only determined by the exact value of 
the valley splitting in real system, which is also the cri- 
terion of the ground state configuration. 

To elucidate the valley and spin configurations of the 
ground state in the absence of the magnetic field, we draw 
a phase-diagram-like picture, Fig.[TJa), by calculating the 
energy difference Ai^J between the lowest level with pure 
("— ") configuration and that with mixed ("m") config- 
uration. From this figure, one finds three regimes for 
each half- well width d, as labeled by (S, Stot)- It is seen 
that the ground state is of "— " ("m") valley configura- 
tion with large (small) valley splitting as expected. For 
the "—"-configuration case, the ground state can vary be- 
tween quartet and doublet as the dot size changes, with 
the transition occurring at around do = 26 nm (shown 
as the vertical dotted line). Since the energy of the low- 
est ("fh") configuration state is always higher than 
that of the lowest "— " ("m") configuration one due to a 
finite valley splitting, the ground state can not be of this 
configuration. However, when A£'q — 0, the lowest "fh"- 
configuration level is degenerate with "m" -configuration 
one, therefore, can also be the ground state in that case. 

In our calculation, we determine the borderline be- 
tween different valley configurations from the degenerate 
condition between them. Specifically, we calculate the 
energy of the lowest "— " and "m" levels in the absence 
of the valley splitting, i.e., and i?™ separately, and ob- 
tain A^;^ = El-E^- AE^. The valley splitting on the 

borderline is then given by AEq = E^ — E^. Actually, 
the real value of the valley splitting A£'q for certain well 
width d should be determined by A J q in Appendix Kl If 
the coordinate (do, A_Eq) locates above the curve of the 
corresponding well width, i.e., AEq > AEq or AE^ < 0, 
the ground state is of "— " valley configuration. Other- 
wise, it should be of "m" -configuration. 

The origination of the variation of the ground state 
spin configuration in the "—"-configuration regime is 
found to be the competition between the orbital and 
Coulomb energies. The relative orbital (Coulomb) en- 
ergy between the lowest quartet and doublet AEq {AEc) 
is defined as the orbital (Coulomb) energy of the lowest 
quartet subtracting that of lowest doublet. One calcu- 
lates the orbital energy from the expectation of Y^i=i ^1 
of the three-electron eigenstates, while the Coulomb en- 
ergy is similarly given by + + . From the 
explicit expressions of the quartet and doublet given in 
Appendix |B1 one notices that all three electrons in quar- 
tet must occupy different orbits while two electrons in 
doublet can stay in the same one. Therefore, the lowest 
quartet states have higher orbital energy than the lowest 
doublet states (A£'o > 0). We finds that the Coulomb in- 



teraction of the lowest quartet states is smaller than that 
of the lowest doublet states {AEc < 0). In Fig.[T](b), we 
plot AEo and —AEc as function of effective diameter, 
where an intersecting between these two quantities can 
be seen. As a qualitative understanding, AEq oc 1/do 
can be characterized from Hujq = h^irHrritd^), while 
AEc oc l/do is estimated from e'^/(47reo'«do)- Therefore, 
the orbital energy is more sensitive to the diameter than 
the Coulomb interaction. For a small diameter QD, the 
relative position between the lowest doublet and quartet 
is dominated by the orbital energy and the ground state 
is doublet. As do increases, AEo can become smaller than 
—AEc with the crossover at do ~ 26 nm in Fig.[T](b), re- 
sulting in the transition of the ground state from doublet 
to quartet. 

From Fig. [TJa), one can directly read the ground state 
configuration with the knowledge of the dot size and val- 
ley splitting, e.g., Sil-Si4 shown as crosses. Moreover, 
one notices that our phase-diagram-like picture is robust 
against the well width, resulting from the strong confine- 
ment regime along z-direction (d ^ dp). 



3. Three- electron spectrum with "— " valley configuration 

In this part, we take advantage of our model to cal- 
culate the three-electron energy spectrum with "— "- 
valley configuration. We study the perpendicular mag- 
netic field dependence at d = 2.12 nm (corresponding to 
A£^Q « 1.9 meV) with do = 20 nm and 29 nm, where the 
ground states lie in the (—,1/2) and (—,3/2) configura- 
tion regimes, respectively. The lowest few energy levels 
are plotted as function of B± in Fig. These eigenstates 

are denoted as (9L3/21 '5i/2' Q3/2' -^-"1/2 ^'^'^ 

-D|y2 1 according to the total spin Stot and Sz of the major 
components in these states. The lowest doublet (quartet) 
states at B± = T are labeled as circle (open triangle). 
It is seen that the ground state for do = 20 nm is four- 
fold degenerate doublet, while that for do = 29 nm is 
quartet, consistent with Fig.[lja). As the magnetic field 
increases, the degeneracy is lifted due to the effect of the 
magnetic field on Landau level and Zeeman splitting, re- 
sulting in many intersecting points as shown in Fig. HJa) 
and (b). We find that some intersecting points show an- 
ticrossing behavior due to the SOCs.® In this work, all 
the anticrossing points are labeled by open squares. In 
the inset, we enlarge the spectrum in the vicinity of one 
anticrossing point at B±^ ~ 0.72 T, where an energy gap 
^ 0.24 /xeV is present ji^Ti^i As reported, the anticrossing 
points are important for spin manipulation, due to the 
strong spin mixed at these points^^ ^^'^^" — 

Since the eigenstates are strongly mixed, e.g., the 
largest component in the ground states is usually in the 
order of 20%, one can not naively describe an eigenstate 
by a single basis function. However, away from the an- 
ticrossing points, the SOCs are weak and the total or- 
bital angular momentum L = Zi -t- ^2 + ^3 is still a good 
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FIG. 2: (Color online) The lowest few energy levels vs. per- 
pendicular magnetic field B±_ in single QDs with d = 2.12 nm 
and the effective diameter do = 20 nm (a) and 29 nm (b). 
The anticrossing points are labeled by open squares. The 
circle and open triangle at B±_ = describe the lowest dou- 
blet and quartet states separately. The total orbital angular 
momentum L = l\ -\- I2 + h of the energy level are labeled 
correspondingly.— The inset in (a) zooms the vicinity of the 
anticrossing point at B±_ ~ 0.72 T. The inset in (b) enlarges 
the spectrum in the vicinity of .B ~ T. 



quantum number. In Fig. [2ja), we labeled this quan- 
tum number for each state. With this quantum number, 
one can understand why an intersecting point is anti- 
crossing or not in the presence of SOCs. By rewriting 
the SOCs in the form of ladder operators, one obtains 
= -ii^^p+s- - P-S+)~^{P+S+ +P-S-) with 
P± = {P^ ± iPy)/2 and ^ S^± iSy. Since P± (5±) 
changes L (Sz) by one unit, the states with (L, Sz) can 
only couple with the states with {L^l, 6*2 ±1) via Rashba 
term and those with {L ± 1, 5*2 ± 1) via IIA term.^ For 
example, in the inset of Fig. Ufa), the spin-down doublet 
with (1, — i) (black chain curve) couples with the quartet 
with (0, ^) (purple dashed curve) via IIA term, resulting 
in the anticrossing behavior. 

To exclude the orbital effect, we apply a parallel mag- 



FIG. 3: (Color online) The lowest few energy levels vs. par- 
allel magnetic field -B|| in single QDs. Open squares indi- 
cate the anticrossing points. The circle and open triangle at 
B± — describe the lowest doublet and quartet in the absence 
of the magnetic field. The total orbital angular momentum 
L — h + I2 + h are labeled. The half- well width d = 2.12 nm 
and the effective diameter do = 20 nm. 

netic field instead. The results with d = 2.12 nm and 
do — 20 nm are shown in Fig. O where the eigen- 

states are denoted as Q 3/2' ^ 1/2' ^1/2' ^3/2' ^ 1/2 

and -D|/2 according to ^tot ^^nd the x-direction com- 
ponent Sx of their major components. In this case, 
the energy levels are linearly dependent on the mag- 
netic field due to the Zeeman splitting. The degener- 
acy from the orbital degree of freedom survives. For ex- 
ample, the ground state for -By < 1 T is two-fold de- 
generate with L = ±1 as shown in Fig.[3l Similar to 
the situation with the perpendicular magnetic field, we 
also observe anticrossing points, e.g., the ones marked 
as open squares at -By ^ 1.1 T. However, we should 
point out that the criterion of the anticrossing points 
here are different from that with perpendicular mag- 
netic field.- In this case, the SOCs can be written as 
Hso = [ao{P+ +P-)- iho{P+ - P")] {S+ + S-)/h with 

S"^ — Sy i %Sz^ 



4. Three- electron spectrum with "m" valley configuration 

As shown in TableUand Fig.[Tl the valley configuration 
of the ground state can be "m" for small valley splitting 
case, e.g., in Si4 with ^ 0.12 /ieV. In this part, we study 
the energy spectrum with the "m" valley configuration to 
illustrate the role of the valley degree of freedom in the 
three-electron spectrum. In Fig. [51 we plot the perpen- 
dicular magnetic field dependence of the energy spectrum 
with the "m" valley configuration in Si4 (d = 3.926 nm 
and do = 30.3 nm). Similarly, the lowest several eigen- 

states here are denoted as Q^Z^j^, Q!~l/2' ^1/2' ^3/2' 
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D 



(-) 



-^^2 -^1/2' according to the total spin Stot and 
of the major components. The lowest doublet and 
quartet states at B± = T are labeled as red circle and 
open triangle, respectively. In the "m" valley configu- 
ration case, each single-electron orbit (distincted by the 
quantum numbers n, I and n^) is four- fold degenerate 
due to the spin and valley degrees of freedom. Therefore, 
at most two electrons in quartet can stay in the same 
single-electron orbit, while all the three electrons in dou- 
blets can stay in the same one. This gives rise to the 
difference between energy spectrum of the "m" and "— " 
valley configurations. 

Interestingly, we find that the magnetic-field indepen- 
dent degeneracy between quartet and doublet can exists 
in the "m" valley configuration. We notice that the de- 
generacy due to the valley degree of freedom has been 
discussed by Wang et al. in two-electron case in Ref.H. 
Here, we find that the present degeneracy can be either 
two-fold or three-fold. In Fig. SI the three-fold degener- 
ate levels (one quartet state and two doublet states) are 
labeled by crosslets, while the two-fold degenerate ones 
(one quartet state and one doublet state) are marked 
by solid triangles. Actually, this energy-level degener- 
acy mainly results from the negligibly small intervalley 
Coulomb interaction and SOCsjii^ One finds that the 
doublet states in the two-fold degenerate case are mainly 
comprised of lU^II,^'') elements. This is associated with 

the situation that 1 15^12^ '') basis functions of the "m" 

configuration are decoupled with |£'2i/2^'') ^^"^ l-^±i/2^^) 
in the absence of the intervalley Coulomb interaction and 
SOCs. Moreover, one finds that there is one \Q^^if2) 
sis function with the same orbital construction as each 



\D 



(m)(2)> 



156 



The 



'±1/2 ) ^'^"^ ^'^ce verse, according to Appendix[ 
corresponding Hamiltonian elements for these two sets 
are equal, i.e.. 



,p,(m)(2)|rr \-ni(m){'2)\_,i^(m) ,„ 

(iViiV2iV3|i/tot|iV{iV^iV^) ± \g^lBB, 
-{NiN2N:i\Htot\N!,N[N'^) 



tot 



(9) 



Thus it's clear that each quartet energy level with Sz = 
±i degenerates with one doublet state (constructed by 

I n(™)(2)\A 

|-^±l/2 I)- 

For the three-fold degenerate levels shown as the curves 
with crosslets in Fig. 01 we find the additional doublet 



state is combined by iD^^^j^^) and l-D^'^^j'*'') basis func- 
tions. However, the combination of these basis functions 
is very complex, therefore, we skip the discussion in this 
work. 



(m)(3)\ 



B. Vertically coupled double quantum dots 

In this part, we investigate the vertically coupled DQD 
case with a perpendicular magnetic field = 0.1 T with 




35.82 



FIG. 4: (Color online) The lowest few energy levels with the 
"m" valley configuration vs. perpendicular magnetic field Bx 
in single QDs. Open squares indicate the anticrossing points. 
The curves with crosslets denote the three-fold degenerate 
levels (two doublet states and one quartet state) and the ones 
with solid triangles represent the two-fold degenerate levels 
(one doublet state and one quartet state). The circle and 
open triangle at B± = describe the lowest doublet and 
quartet and doublet. The total orbital angular momentums 
L = l-i-\-l2 + l3 are labeled. The half-well width d — 3.926 nm 
and the effective diameter do = 30.3 nm. 



both the "— " and "m" valley configurations. We first 
study the interdot barrier-width dependence of the "— " 
valley configuration. In Fig. [Sfa), the orbital energies of 
the lowest two subbands along the z-direction are plotted 
as function of interdot barrier width. One can see that 
the orbital energy of the first subband increases while 
that of the second subband decreases as the barrier width 
increases, resulting the decrease of the energy difference 
El — Eq between them. This corresponds the decrease of 
the interdot coupling. Differently, the valley energies of 
the lowest two subbands present oscillation against the 
barrier width, which can be expected from the formula 
given in Appendix [Xj Since the intersubband coupling 
due to the valley-orbit coupling is much weaker than the 
intersubband splitting from the orbital degree of freedom, 
we neglect this term in our calculation. Actually, the first 
subband dominates the lowest few states because of the 
large intersubband splitting (> 10 meV) in the strong 
interdot coupling regime studied here. The lowest four 
levels are plotted in Fig.[5Jc), where the energies from 
the subband and valley are subtracted for improving the 
solution of levels. Similar to the single QD case, the 
quartet-doublet transition of the ground state can occur 
and the anticrossing behavior (illustrated as open squares 
at a ^ 1.5 nm) can be observed. We also study the 
barrier-height dependence with "— " valley configuration 
and plot the results as Fig. [Sl[e)-(g). In Fig. [Sfg), the 
quartet-doublet transition of the ground state and the 
anticrossing properties are also present. 
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FIG. 5: (Color online) (a) and (e) represent the interdot 
barrier-width and barrier-height dependences of the lowest 
two subband energies along z-direction. (b) and (f) are the 
valley energies of the lower valley eigenvalue of the lowest two 
subbands as function of the interdot barrier width and barrier 
height, (c) and (g) show the three-electron energy spectrum 
with "— " valley configuration, while (d) and (h) give the cor- 
responding spectrum with "m" valley configuration. Open 
squares in (c) and (g) indicate the anticrossing points. The 
curves with crosslets denote the three-fold degeneracy of the 
levels (two doublet states and one quartet state). The total 
orbital angular momentum L — h + I2 + I3 of each states are 
labeled. Here, the effective diameter do = 29 nm, the well 
width d = 2 nm and Bx = 0.1 T. In (a)-(d), the barrier 
height Vo = 50 meV. In (e)-(h), the barrier width a = 1 nm. 



IV. SUMMARY 



In summary, we have studied three-electron energy 
spectra in Si/SiGe single QDs and vertically coupled 
DQDs by using the exact-diagonalization method. In 
our calculation, the Zeeman splitting, SOC, valley cou- 
pling and electron-electron Coulomb interaction are ex- 
plicitly included. Due to the strong Coulomb interaction, 
a large number of basis functions are employed to con- 
verge the energy spectrum. The ground state energies in 
single QDs show good agreement with the experiment. 
As a supplement of the experimental data, we identify 
the valley configuration, spin configuration and the de- 
generacy factor of the ground state from our calculation. 
We then systematically study the ground state configu- 
ration in the absence of magnetic field, and find that the 
ground state is of pure and mixed valley configurations 
with large and small valley splittings, respectively. We 
show that the ground state with mixed valley configu- 
ration is always doublet. In contrast, the ground state 
with pure valley configuration is doublet in small dots, 
while it can be quartet in large dots. Then, we explicitly 
study the energy spectra of the three-electron states with 
two typical valley configurations, i.e., pure and mixed 
valley states, in the presence of perpendicular and par- 
allel magnetic fields. In the pure valley configuration 
case, we show that the doublet-quartet transition of the 
ground state can be realized by tuning the magnetic field 
and/or dot size. For the mixed valley configuration, the 
doublet-quartet transition of the ground state can also be 
realized by magnetic field and the three-electron energy 
spectrum presents interesting quartet-doublet degener- 
acy, which is connected with the negligible intervalley 
coupling. Due to the spin-orbit coupling, the intersect- 
ing points between the energy levels with identical valley 
configuration but different spin quantum numbers can 
present anticrossing behavior, which is expected to bene- 
fit the manipulation of spin. We point out and analyze all 
these anticrossing points in each case. Furthermore, we 
study the barrier-width and barrier-height dependences 
of the three-electron energy spectrum in vertically strong- 
coupled DQDs. We also find anticrossing points in the 
pure valley configuration case and quartet-doublet de- 
generacy in the mixed configuration case. 
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Appendix A: VALLEY COUPLING ELEMENTS 

^ 1^1^VGo! + 2^.^06-0^ sm2(fcod), (Al) 
.1 _ K^"fco'Co'cos[2fcsi(f 

+ 2V;K)Co^ sin2(fcofi) cos(fcsia), (A2) 

A?! = M!^l!^ + 2K^oCi2sin2(fcid), (A3) 

1 _ K^'fci'Ci'cos(2fcsi(f 
^1,1 — 

+ 2KVbCi2 sin2(fcid) cos(fcsia), (A4) 

Al Al . r Kfi^^O^lCoCl ■ r, / , oJM 

= -Aio = i{ sinfcsia + 2(i) 

+ 2K^^oCiCo sin(/cod) sin(/cid) sin(/csia)}, (A5) 

where Vv = 7.2 x 10^^^ m represents the ratio of the 
vahey couphng strength to the depth of the quantum 
weU potential^ 

Appendix B: THREE-ELECTRON BASIS 
FUNCTIONS 

The three-electron basis functions can be constructed 
in form of either doublet {Stot = 1/2) or quartet (5tot — 
3/2). We first combine the spin wavefunctions via 
Clebsch-Gordan coefficients. For the doublet states, the 
spin functions can be expressed by 

Xi/2 = [(^"2 + ai/32)Q;3 - 2aiQ;2/33]/V6 

x\/2 = [2/3i/32a3-(/3ia2+ai/32)/33]/V6 

'Xa/2 = (/3ia2 - ai^2)a3/V2 

x'^1/2 = (/3ia2 - ai/32)/33/V2, (Bl) 

and for the quartet states, one has 

xf/2 = aia2a3 

xf/2 = (aia2/33 + ai/32a3 + /3ia2a3)/V3 
X-1/2 = ("1/32/33 + /3i/32a3 + /3ia2/33)/\/3 
x!3/2 = PihPz. (B2) 
I 



Here, ai = \ t)i and jSi = \ with the subscript de- 
noting the three electrons. The superscript "S" labels 
the symmetric configuration for the permutation of any 
two electrons, while "A" and "p" represent the symmet- 
ric and antisymmetric configurations for the permutation 
only for the electrons "1" and "2", respectively. With 
these spin wavefunctions, we then add the correspond- 
ing orbital parts to construct the total basis functions. 
When the three electrons occupy three different orbits, 
i.e., Ni 7^ Nj with z, j e {1, 2, 3}, one has the total wave- 
functions of the doublet states 

- 0'"^v^/2 - 0^<^^x?/2. 

m^^) - 0''<^)x^/2-0^(^'x^/2, 
\D^^lf^) = 0''(^)x^/2-'^'^^^x^/2, (B3) 

and those of the quartet states 

IQ3/2) = '/'''xf/2, 

\Qf,{) = '/'^xf/2, 

|Q%> = <l^\U/2- (B4) 

When two electrons with opposite spins share the same 
orbit, i.e., Ni = N2 ^ N3 or Ni ^ N2 = N3 or iVi = 

N3 ^ N2, 

l^l/f ) = 0''xt/2-0'x?/2, 

l^i-i^;2^) = 0''x^i/2-'/''x^_i/2- (B5) 

Here, {(f)} are the corresponding orbital wavefunctions. 
The superscript "A" on the right hand side of Eq. (jB4p 
describes the permutation antisymmetric character for 
the exchange of any two electrons and the "(1)" and "(2)" 
in the superscripts of Eq. (jB3|) are used to distinguish the 
two doublet configurations with the same spin magnetic 
quantum numbers or the two orbital functions with the 
same symmetry. The orbital wavefunctions can be ex- 
pressed as 



</>^ - -1={\N,N2N,) ^ 
^{\N,N2N,) 
^{\NiNsN2) 
^(2\N^N2Ns] 

^(l^3A^liV2) 



0P(i) 
6A(i) 

6P(2) 



\N1N3N2) + \N2N3N1) - \N2N1N3) + liVgiViiVs) - \N3N2N1)) 

- \N2N3Nl) + \N3NlN2) - \NlN3N2)) 

h |^2iV3iVi) + \N3N2N1) + \N3N1N2) - 2\NiN2N3) - 2\N2NiN3)) 

- 2\N2NiN3) - \N2N3N1) + \N3N2N1) - \N3N1N2) + \N1N3N2)) 

- \N3N2N1) + \N1N3N2) - \N2N3N,)), 



(B6) 
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i(|7ViiV3iVi) 

i(|7ViiV2iV2) 

U\NiN2Ni) 



\N3N1N1)), 
\N2N,N2)), 
\N2N,N,)), 



Ni=N2^ N3, 
Ni^N2 = N3, 
Ni^Ns^ N2, 



(B7) 



^{2\N,N,N,) 
^{\N2N^N2}^ 
^,miNiN2} 



\N1N2N2) - 2\N2N2N,)),N, ^N2^ N3, 
- \N1N2N1} ~ \N2NiNi)),Ni =N3^ N2. 



(B8) 



Appendix C: MATRIX ELEMENTS OF COULOMB INTERACTION 

The Coulomb interaction components lead^^^^^ 
{N,N2N,\Hi^\N[N!,N!,) = --^5^^,^,<5,^+,^,,+,, ^ ^n:.'7nl.4,4.^(C;.n„' <^</>L.' 

(CI) 

with {i, J, fc} = {1, 2, 3}. Here, \Ni) = Inilin^invi) is the single-electron wavefunction of the i-th electron. Superscripts 
7j(j) and ^y'^^j-^ run over the two valleys with — 1 and 77^ = —77! = 1 (Ref.181). U in Eq. (jCip can be expressed as^^ 

J J —00 

(C2) 

P^i^ can be obtained from the lateral part of the integral {nlnz'y\e^^''^\n'Vn'^Y) . 

I j ^"j ^2 n' n J 

with ; = ^r- n = i + .7 + and sgn(a:) being the sign function. Wn':^' in Eq. jCU is the integral 
(n^7|exp(jfc^z)|n^7'). ^ 
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From the explicit expression of the three-electron Hamilto- 
nian, one can demonstrate that the Hamiltonian without 
spin-orbit term is non-vanishing only between states with 
the same total orbital angular momentum L = h + I2 + h- 
And the spin-orbit coupling is very weak compared with 
the orbital energy and Coulomb interaction. Therefore, one 
can denote each state (except in the vicinity of the anti- 
crossing points) by the total orbital angular momentum. 
^® Here we take n„i = n„2 = — and n„3 = -I- for the permuta- 
tion antisymmetric or symmetric characters of the orbital 
part of the basis functions for the electrons "1" and "2". 
Even though, the identity among the three electrons is still 
guaranteed by the orthogonality and completeness of the 
entire Hilbert space. 



